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We investigate collective excitations coupled with baryon density in a system of massless three- 
flavor quarks in the collisionless regime. By using the Nambu-Jona-Lasinio (NJL) model in the 
mean-field approximation, we field-theoretically derive the spectra both for the normal and color- 
flavor locked (CFL) superfluid phases at zero temperature. In the normal phase, we obtain usual 
zero sound as a low-lying collective mode in the particle-hole (vector) channel. In the CFL phase, 
the nature of collective excitations varies in a way dependent on whether the excitation energy, ui, is 
larger or smaller than the threshold given by twice the pairing gap A, at which pair excitations with 
nonzero total momentum become allowed to break up into two quasiparticles. For ui <C 2A, a phonon 
corresponding to fluctuations in the 17(1) phase of A appears as a sharp peak in the particle-particle 
( "H" ) channel. We reproduce the property known from low energy effective theories that this mode 
propagates at a velocity of vh = l/\/3 in the low momentum regime; the decay constant fu obtained 
in the NJL model is identical with the QCD result obtained in the mean-field approximation. We also 
find that as the momentum of the phonon increases, the excitation energy goes up and asymptotically 
approaches uj — 2 A. Above the threshold for pair excitations (to > 2A), zero sound manifests itself 
in the vector channel. By locating the zero sound pole of the vector propagator in the complex 
energy plane we investigate the attenuation and energy dispersion relation of zero sound. In the 
long wavelength limit, the phonon mode, the only low- lying excitation, has its spectral weight in 
the H channel alone, while the spectral function vanishes in the vector channel. This is due to 
nontrivial mixing between the H and vector channels in the superfluid medium. We finally extend 
our study to the case of nonzero temperature. We demonstrate how Landau damping smears the 
phonon peak in the finite temperature spectral function. We find a pure imaginary pole of the H 
propagator in the complex energy plane, which can be identified as a diffusive mode responsible for 
the Landau damping. From the pole position we derive the thermal diffusion constant. 

PACS numbers: 12.38.-t, 25.75.Nq 



I. INTRODUCTION 

Quark matter at high baryon density is considered to 
have a rich phase structure in a way dependent on the 
physical conditions for color and flavor such as weak equi- 
librium and color neutrality. For the past few decades 
it has been predicted that cold quark matter is color 
superconducting when the baryon density is sufficiently 
high [lj . Color superconductivity occurs in quark matter 
in a similar manner to ordinary superconductivity in met- 
als; any attractive force between quarks can bring about 
an instability of the Fermi surface against Cooper pair- 
ing, which subsequently leads to a BCS superfluid state 
characterized by condensation of quark Cooper pairs and 
to a color version of the Meissner effect through nonva- 
nishing color charge carried by the pairs. In fact, the 
quark-quark force induced by one-gluon exchange is at- 
tractive in the color antisymmetric channel. At asymp- 
totically high densities where the interactions are domi- 
nated by one-gluon exchange, this attractive force should 
drive Cooper instability. More recently, effective model 
analyses, rather than rigorous QCD calculations, have 
predicted that quark matter is a color superconductor 
even at several times normal nuclear density and that 
the zero temperature pairing gap and thus the critical 
temperature are of order 10-100 MeV. 

In possible physical realizations in heavy ion collisions 



and neutron stars, color superconducting quark matter 
might leave evidences for color superconductivity via pre- 
cursory phenomena above the critical temperature 0, 
neutrino emission from proto-neutron stars 0], and neu- 
tron star cooling 0,0, etc. To obtain such possible evi- 
dences, it is indispensable to clarify the equilibrium prop- 
erties of color superconductors such as the equation of 
state and elementary excitations. This paper elaborates 
on collective excitations in quark matter, which might 
help us detect the onset of color superconductivity. 

Since quarks have more internal degrees of freedom 
than electrons by two (color and flavor), various pairing 
patterns are possible among quarks, even within a state 
of S-w&ve, spin-singlet pairing. The ratio between the 
strange quark mass and the quark chemical potential is 
known to be a crucial parameter that determines the en- 
ergetically favorable pairing pattern. In the case of three 
massless flavors the color-flavor locked (CFL) phase 
is the ground state 0. In this paper we shall focus only 
on the massless limit of CFL quark matter in which the 
baryon number conservation is spontaneously broken by 
the condensation of quark pairs. This can be mathemat- 
ically formulated as global U(1)b symmetry breaking, as 
we briefly recapitulate in Sec. ITT1 

Superfluid 3 He, a typical laboratory fermion system 
that exhibits global U(l) symmetry breaking 0, is sim- 
ilar to CFL quark matter in many respects. In fact, the 
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ground state of liquid 3 He, i.e., the B phase, has a con- 
densate isotropic in spin-orbit space, which is similar to 
the CFL phase in which the condensate is isotropic in 
color-flavor space. In both systems, responses to rota- 
tion are characterized by a lattice of vortices |9| . 

In the B phase of superfluid 3 He in which the pair- 
ing gap has no nodes, thermal quasiparticle excitations 
are suppressed at temperatures much lower than the gap. 
In this case, the thermodynamic properties such as the 
specific heat are presumably controlled by phonon exci- 
tations, Nambu-Goldstone bosons corresponding to fluc- 
tuations in the U(l) phase of the pairing gap, which ap- 
pear in the low energy spectral function in the particle- 
particle channel. Then, the same kind of phonons in CFL 
matter (often denoted by H) are expected to play a key 
role in the thermodynamic properties at low tempera- 
tures 0, [n| . In this paper we investigate the H phonon 
spectrum in the energy range including the scale beyond 
the pairing gap A in the random phase approximation of 
a Nambu-Jona-Lasinio (NJL) model for CFL quark mat- 
ter. This microscopic approach enables us to see that H, 
if embedded in the quasiparticle continuum lying above 
2A, would decay into two quasiparticle quarks. 

We remark that an H mode does not appear as a trans- 
verse mode, a feature responsible for a nonzero value of 
the superfluid baryon density at zero temperature [l2| . 
Transverse excitations will be ignored in the present 
study, but play a role in breaking up phase coherence in 
the superfluid state and hence are relevant for responses 
to magnetic field and rotation Qj, fl3| and for reduction 
in the transition temperature [14j. 

In normal quark matter we also describe zero sound as 
fluctuations in the baryon density by analyzing the spec- 
tral function in the particle-hole (vector) channel in the 
random phase approximation in which collisions between 
quark quasiparticles are ignored (see Ref. 15] for general 
arguments of relativistic zero sound) . We extend the de- 
scription of zero sound to the case of the CFL phase. We 
find that in the vector channel, zero sound manifests it- 
self as an attenuated mode in the quasiparticle continuum 
lying above twice the gap, while H dominates the spec- 
trum below the continuum. We remark that for general 
excitation momenta the random phase approximation is 
strictly applicable to weakly coupled systems. In this 
respect, a recent observation of zero sound in the lattice 
simulation based on a QCD-like model [lf| is noteworthy. 

Note that in the collisionless regime the restoring force 
of the collective modes is a self-consistent field formed 
by a number of particles, which would be disrupted by 
interparticle collisions if any. Collective modes in the 
hydrodynamic regime such as first sound, of which the 
restoring force is provided by the collisions, will be ad- 
dressed in our future study. For the study of low-lying 
collective modes in both regimes, the Landau theory of 
Fermi liquids is useful. This is the case with normal 
liquid 3 He 0, superfluid 3 He 0, and normal quark 
matter [l^. 

In calculating the spectral functions in the CFL phase, 



it is important to take into account the mixing between 
the particle-particle and particle-hole channels. This is 
because such mixing plays a role in determining the spec- 
tral weight of H and zero sound modes. We find that in 
the long wavelength limit the spectral function in the 
particle- hole channel vanishes due to the mixing effect. 

We summarize the main conclusions of this paper as 
follows: 

• The velocity of the H mode in the small momen- 
tum regime is l/v3 at zero temperature. As the 
momentum goes up, the mode energy monotoni- 
cally approaches 2A, above which individual pair 
excitations are kinematically allowed. The velocity 
l/v^ corresponds to the ratio of the spatial and 
temporal components of the decay constant fjj ■ We 
find that the expression for fn in the NJL model 
adopted here is the same as the QCD result in the 
mean-field approximation. 

• In the CFL phase zero sound appears as an attenu- 
ated mode in the quasiparticle continuum extend- 
ing above 2A. The attenuation width decreases 
with increasing momentum. 

• In the CFL phase, the spectral function in the vec- 
tor channel is strongly suppressed at small mo- 
menta by the mixing with the particle-particle 
channel. 

• At finite temperature at which thermally excited 
quasiparticles are present, the H mode has a finite 
decay width due to Landau damping. Correspond- 
ingly, a pure imaginary pole appears in the H prop- 
agator. We estimate the thermal diffusion constant 
from the pole behavior in the complex energy plane. 

We note that the present work, which focuses on col- 
lective modes coupled with baryon density fluctuations in 
the CFL phase,_is a complement to the work by Gusynin 
and Shovkovy 20] that gave a detailed account of collec- 
tive modes coupled with color current fluctuations in the 
CFL phase. 

This paper is organized as follows: We briefly review 
the CFL Cooper pairing and the symmetry breaking pat- 
tern in Sec. [HJ We describe a model for quark matter, 
the self-consistency equations, and the spectral functions 
in Sec. IIIII Full expressions for the spectral functions are 
listed in Appendix lAl Section Hvl is devoted to showing 
the results for H and zero sound modes at zero temper- 
ature. In Sec. EJ we present the results at finite tem- 
perature and discuss Landau damping of the H mode. 
Concluding remarks are given in Sec. IVII Appendix [B] 
explains how we perform the analytic continuation of 
the propagators of collective excitations to search for the 
poles in the complex energy plane. We work in units 
ft, = c = fee = 1- 




II. COOPER PAIRING, SYMMETRY 
BREAKING, AND NAMBU-GOLDSTONE 
BOSONS 

In this section, we summarize the fundamental prop- 
erties of the CFL phase. In particular, the symmetry 
breaking pattern and the associated Nambu-Goldstonc 
bosons are mentioned. 

In the CFL phase with itds-fiavor and RGB-co\oy 
massless quarks, the energetically favored Cooper pairing 
of quarks is antisymmetric in color space and has even 
parity and zero total angular momentum |y, |jj • Then, 
the antisymmetry of the pairing in flavor space follows 
from the Pauli principle. [Note that in the absence of 
quark masses, only quarks of the same chirality can be 
paired |2l| .] The corresponding diquark condensate reads 

(^75^) ~ A {S m S Pj - SajSpi), (1) 

where the Roman and Greek subscripts stand for the 
indices in flavor and color space, respectively, ip is the 
quark spinor, and ip c = CifF is the charge conjugate 
spinor with the charge conjugation matrix C — i^j 
in the Pauli-Dirac representation. A precise relation be- 
tween the condensate and the pairing gap will be given in 
Eq. (J5J after introducing a model for quark matter. The 
coexisting condensate in the CFL phase, which arises in 
a color and flavor symmetric state, will be ignored since 
it makes negligible difference in the present analysis that 
focuses on the superfluid properties. 

The condensate Q breaks both SU(3) l and SU (3)_r of 
chiral symmetry, whereas it leaves unbroken a symmetry 
under a vector rotation in flavor space and a simultaneous 
color rotation in the opposite orientation. The symmetry 
breaking pattern in the CFL phase can thus be written 
as HH 

[SU(3) color ] x SU(S) L x SU(3) R x U{1) B 
- SU(3) C+L+R x Z 2 . 

Here, Z 2 is a symmetry under ip — > —ip. Nine Nambu- 
Goldstone bosons result from the spontaneous symmetry 
breaking described above. Eight pion-like bosons associ- 
ated with the breakdown of the axial part of chiral sym- 
metry form an octet in SU(3)c+l+r just like (ir,K,rj) 
in hadronic matter. One more massless boson arises as 
a color-flavor singlet from the spontaneous breaking of 
baryon number conservation in a manner that preserves 
local electric and color charge neutrality. This singlet bo- 
son is a phonon (often referred to as H) corresponding 
to U(1)b phase fluctuations of the condensate Q). 

For completeness, let us consider the other phases 
which are predicted to appear in the phase diagram 
for neutral quark matter with nonzero strange quark 
mass [23.1231124) and confirm that all these phases do not 
have an H mode. It is well known that in the 2SC phase 
in which only two colors and two flavors participate in 
pairing, there remains an unbroken baryon number sym- 
metry associated with the unpaired color and flavor. In 



FIG. f: Schematic picture of the H propagation. Colored 
diquark condensates supply color charge in such a way as to 
make H a color singlet at the edges. 

the same way, the uSC and dSC phases in which only ds 
and us pairing is breached as compared with the pair- 
ings in the CFL phase, respectively, preserve a modified 
U(1)b symmetry. The uSC phase, for instance, has fi- 
nite condensates (ud) and (su). These two condensates 
are invariant simultaneously under rotations generated 
by an appropriate linear combination of the U(1)b gen- 
erator and two U(l) generators corresponding to flavor 
number conservation. Therefore, there is no H in the 
2SC, uSC, and dSC phases. Put another way, in each 
of these phases, U{\)b phase fluctuations in the order 
parameter are inevitably coupled with electric and color 
charge density since the condensates, as a whole, are not 
charge neutral in electricity and color. We note that in 
the CFL phase modified by quark masses different among 
flavors, H modes could be affected by coupling with elec- 
tric and color charge density. In this case, the quark 
system, if being color neutral and being neutralized and 
/3-equilibrated by an electron gas, would generally have a 
superfluid part that is no longer charge neutral in electric- 
ity and/or color. This nonneutral superfluid part would 
occur not only at nonzero temperature, but also in the 
presence of gapless quark modes even at zero tempera- 
ture m. 

All the Nambu-Goldstone bosons in the CFL conden- 
sate {Q} are color singlets because the corresponding fluc- 
tuations in the order parameter are not coupled with 
color charge density and thus are not eaten by the lon- 
gitudinal component of color gauge fields through the 
Anderson-Higgs mechanism. It follows that the CFL pi- 
ons and H consist of at least four and six quarks, re- 
spectively j^Qj, i.e., they are composed of such combina- 
tions of particle-particle and hole- hole pairs as (RG) (RG) 
and (RG)(GB)(BR). In this paper we assume that col- 
ored diquark condensates compensate color charge at the 
edges of the H propagator as shown in Fig. We can 
then describe the propagating part of H by a state hav- 
ing an excited quark Cooper pair in Fock space. This 
picture is based on a mean-field approximation, which is 
valid as long as the size of the mean-field A is larger than 
the size of quantum and thermal fluctuations around A. 



III. MODEL CALCULATIONS 

In this section we first describe the equilibrium prop- 
erties of a normal fluid and a superfluid of quark mat- 
ter with three massless flavors by utilizing an effective 
model with local four-quark interactions. We then ex- 
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amine excitation spectra in both phases in the random 
phase approximation. 

A. Self-consistency equations 

We begin with a model for massless three-flavor quark 
matter. For the later purpose of describing collective 
excitations in the particle-particle (H) and particle-hole 
(vector) channels, it is convenient to deal with the La- 
grangian density having the diquark and vector interac- 
tions. One can generally obtain such interaction terms 
from any local four-quark interactions after appropriate 
Fierz transformation. Then the Lagrangian density can 
be written as 

C = $(i$ + nrfty + CD + Cv, (2) 

where fi q is the quark chemical potential, Cd and Cy 
are the diquark and vector interaction terms as will be 
specified below. As we will see later, /j, q receives a finite 
correction from the vector mean-field. Hereafter we will 
take n q — 500 MeV. The corresponding baryon density 
is, as we will see later, approximately lOpo in our calcu- 
lation where po ~ 0.16 fm -3 is the normal nuclear den- 
sity, and may be relevant in the cores of compact stellar 
objects, although ideally symmetric CFL quark matter 
as considered here is only relevant at much higher den- 
sities where the strange quark mass is negligibly small 
compared with the quark chemical potential. Note, how- 
ever, that we have technical difficulty in taking an ex- 
tremely large value of because the effective model 
with the interaction terms that will be shown below is 
non-renormalizable and needs a finite cut-off parameter, 
A. We shall choose A = 1 GeV in this paper. Then, our 
choice, fx q = 500 MeV, can be considered effectively as 
the upper limit at which quark excitations around the 
Fermi surface would hardly suffer from any cut-off arti- 
fact. We thus expect that our analysis of collective exci- 
tations with momenta up to ~ 100 MeV is free from such 
artifact. 

In this paper, for simplicity, we choose Cd to be the 
interaction that occurs only in the color-flavor singlet di- 
quark channel. This ansatz for Cd makes no difference 
in the resulting gap equations as far as ideally symmetric 
CFL quark matter is concerned. The four-quark interac- 
tions thus take the following forms: 

C D = ^*iP%^)(^ H <P^'f), (3) 
C v = -^(^)(^), (4) 

with = vy^e i3 e a pk m which the sum is taken over k. 

The value of Gy is of the order of G, but our quali- 
tative results for elementary excitations are independent 
of a particular choice of Gy. Hereafter we shall take 
Gy = G/2 for definiteness. A different Gy provides 
a different (dimensionless) spectral weight in the vector 



channel, which is approximately scaled by Gy unless Gy 
is anomalously large. 

It should be emphasized that the vector interaction 
generally arises since chiral symmetry is intact; Cy is 
chirally symmetric itself. For specific values of Gy the 
phase diagram on the temperature versus chemical po- 
tential plane can have such intricate structures as to con- 
tain multiple critical end-points associated with diquark 
and chiral condensates [23. Since the density region of 
interest here is presumably beyond the chiral transition 
density, we simply ignore effects of chiral condensation. 

As to the diquark interaction, we shall fix G in such 
a way that the solution to the self-consistency equations 
with Gy = G/2 yields A = 25 MeV with /z, = 500 MeV 
at zero temperature. As we will see in Sec. lIVI the prop- 
erties of collective excitations in the CFL state change 
drastically according to whether they reside above or be- 
low 2A, which corresponds to a threshold above which 
pair excitations are allowed to occur. As long as A 
is much smaller than the quark chemical potential and 
much larger than the temperature, the energy and mo- 
mentum of collective excitations are essentially scaled by 
A. In this case, a single example presented here with the 
choice A = 25 MeV suffices for us to deduce what takes 
place for general values of A. 

Let us proceed to adopt the mean-field approximation 
by introducing variational variables, i.e., the gap param- 
eter, A, and the quark number density, Tlq j cLS 

A = f<^0/^>, (5) 

n q = $ 7 ». (6) 

As we shall see, we obtain the self-consistency equations 
that determine these variables by substituting the mean- 
field quark propagator with A and n q into the right side 
of Eqs. (J5J and @). Equivalently, these self-consistency 
equations can be obtained from the stationary conditions 
on the mean-field thermodynamic potential with respect 
to A and n q . The corresponding equations are the gap 
equation and the relation between the quark density and 
chemical potential. 

It is well known that the mean-field quark propagator 
in the CFL phase takes the simplest form in the CFL 
basis, ip ai = (\ A /V2) ai ii A and 4> al = xp A (X A /\/2) ia . 
Here, the matrices, X A , in color-flavor space are th e Ge ll- 
Mann matrices for A = 1, . . . , 8 and (\°)i a = y/2/35i a 
for A = 0, which are normalized as ^tr\ A X B = S AB . In 
the CFL basis, the gap matrix in the quark propagator, 
AP'l, has only diagonal components, i.e., i7sdiag{A^}, 
where A^ = — A for an octet in color-flavor space (^4 = 
1, . . . , 8) and A^ = 2A for a singlet (A = 0). 

In the present analysis of the self-consistency equa- 
tions, one is not allowed to simplify calculations by drop- 
ping the antiparticle part in the quark propagator. This 
is because the cancellation between the particle and an- 
tiparticle contributions for momenta above the Fermi sur- 
face is crucial for the evaluation of n q in the NJL model. 
At zero temperature, as a first approximation, one can 
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estimate n q from the number density in a free quark gas 
as 3^/tt 2 = 4.9 fnT 3 for fi q = 500 MeV, which corre- 
sponds to a baryon density of roughly 10po as we have 
mentioned before. The NJL model calculation including 
both the particle and antiparticle contributions results 
in n q = 5.1 fm~ 3 for fi q = 500 MeV with A = 25 MeV. 
This is consistent with the first estimate, while we ob- 
tain n q = — 14.7 fm -3 in the absence of the antiparticle 
contribution; even the sign is inconsistent. 

It is a usual technique to decompose the quark prop- 
agator into the particle and antiparticle parts by using 
the energy projection operators of noninteracting mass- 
less quarks, = ^(1 ± 7°7 • p), where p — p/p. Using 
these operators, we can explicitly write down the mean- 
field quark propagator in 2 x 2 Nambu-Gor'kov space as 



odd n. The results read 



S AB (p») = *ab7° 
where the particle part is 



A- 
A+ 



S A - 



A+ 
o Ar 



(7) 



the antiparticle part is 



P0 + (P-Hr) -i&Al5l° 



(8) 



Po - (P+Mr) -«A A 7 5 7° 
-«Aa757° Po + (p+Mr) 



and we have chosen the Nambu-Gor'kov basis as 



(9) 



4> c {-p») 



(10) 

Here we have defined the energy of quark and antiquark 
quasiparticles as e^ A = (p± /i r ) 2 + A\ and the quark 
chemical potential renormalized by the vector interaction 
term as 



— tanhf ^- ) + — tanhf ^ 



2T / 



C 2A 



n «-^2 I d PP 2 



+ (same with e 




8 f A tanh^V 




dfi r \2TJ 


dfl r 


+ (same with e~ 





2T 
0, (12) 



tanh 



Z 2A 

2T 



= 0. (13) 



The first term in the square brackets is the particle contri- 
bution from octet quarks; the second from singlet quarks. 
The antiparticle contributions lead to the same expres- 
sions as the particle contributions with e~ replaced by 
e+. 

We determine G and fi r from Eqs. fTTJl . lfT2*|) . and ltT3|) 

in such a way that the solutions to these equations yield 
A = 25 MeV with ji q = 500 MeV at zero temperature, as 
we have explained above. Once we fix the model param- 
eters at zero temperature, we can derive the temperature 
dependence of A and \i r from the same equations. Such 
self-consistency equations predict a second-order phase 
transition from the CFL phase to unpaired quark mat- 
ter at T c = 18.15 MeV. This value of T c is larger than 
the BCS value of ~ 0.57A due to the two-gap structure 
of the quark excitations (—A and 2A for octet and sin- 
glet quarks, respectively) j2^|. We note, however, that 
in a situation in which CFL quark matter behaves as 
a type-I color superconductor, thermal fluctuations in 
the gauge fields could play a role in changing the phase 
transition from second to first order and in lowering the 
transition temperature The analysis of this situa- 

tion is beyond the scope of the present paper. In Sec.lVl 
we will present the finite temperature results showing 
an appreciable Landau damping effect on H modes at 
T = 0.8T C , under the expectation that CFL quark mat- 
ter at T = 0.8T C is not affected severely by thermally 
fluctuating gauge fields. 



B. Collective excitations 



fl r = fl q - Gv ■ n q , (11) 



which amounts to fi r = 434.7 MeV for our parameter 
choice at zero temperature. 

We can finally obtain the self-consistency equations at 
given temperature T by substituting the off-diagonal and 
diagonal components of Sab given by Eq. J7J into the 
right side of Eqs. and @, respectively. Here we set 
Po = inirT and take the Matsubara frequency sum over 



We proceed to describe longitudinal collective excita- 
tions in a superfluid as well as in a normal fluid of quark 
matter in the random phase approximation. The col- 
lective modes that we can describe in this approximation 
are an H phonon in the particle-particle channel and zero 
sound in the particle-hole channel, and they lie in the col- 
lisionless regime. 

In the random phase approximation, as we have briefly 
noted in Sec. [HJ we can assume that for H modes, color 
singletness is achieved at the edges of the propagator (see 
Fig. ^| . We can thus regard such collective excitations 
as an excited Cooper pair with nonzero total momen- 
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express the bubble diagrams as 



FIG. 2: Schematic diagrams representing the propagation of 
collective excitations. 



turn that multiply scatter with each other in a super- 
fluid medium, instead of considering the propagation of 
three excited Cooper pairs. We can likewise consider zero 
sound in terms of the propagation of a quark particle-hole 
pair both in a normal and a superfluid medium. Conse- 
quently, the propagator of collective excitations can be 
expressed as a sum over the bubble diagrams in the cor- 
responding channel as exhibited in Fig. |3 In this subsec- 
tion we construct the vertices from the interaction terms 
in the H and vector channels and compute the bubble 
diagrams. 

Let us now construct matrices, T, at the vertices from 
the interaction terms. In doing so, we first rewrite the 
four-quark interactions in terms of the Nambu-Gor'kov 
basis. This is straightforward although we should pay 
attention to how to count the Nambu-Gor'kov replicas. 
Since collective excitations couple not only with the orig- 
inal fields but also with the replica fields, the four-quark 
interactions in terms of the Nambu-Gor'kov basis involve 
the original ones, the replica copies, and the original- 
replica cross terms. The resultant diquark interaction 
is 



C D 4 x C D 
G 
~ 12 



where 



y+ij . 
1 Da/3 ' 



pij n 



p— *J _ 
1 Da/3 ~ 



iP% 



(14) 



(15) 



are the matrices associated with fluctuations in the am- 
plitude and the U(l) phase of the diquark condensate, 
respectively. We can likewise express the vector interac- 
tion as 

£ y ^4x/; v = -^($r^)(*r^) (16) 



with the matrix defined by 



pM 

1 v 



y 
-7" 



(17) 



By using the quark propagator JJJ) with the defini- 
tions iJHJ and © and the matrices (|15|l and i|17|) . we can 



i f T d 4 P 



2 J (2tt) 



4 tTT D S(po+qo,p+q)T D S(p ,p) , (18) 



n y I/ (<7o,<7) 



1 f T d 4 p 

2 J (27T) 1 trT v S (Po + qo,P+q)TvS(p ,p) , (19) 



n M0?o,9) 



tTT^S(p +qo,P+q)T^S(p ,p) (20) 



2 J (2tt) 4 d 

for the H, vector, and -ff-vector mixed channels, respec- 
tively. Here, J T c? 4 p/(27r) 4 denotes the integration over 
p with the cut-off A and the Matsubara summation with 
respect to po, and the trace is taken over color, flavor, 
spinor, and Nambu-Gor'kov indices. The factor 1/2 is 
necessary for adjusting the duplicate loop counting in 
the Nambu-Gor'kov basis. 

In the present study we do not consider the bubble 
diagrams including the vertex associated with fluc- 
tuations in the amplitude of the order parameter. This 
is partly because collective amplitude modes are embed- 
ded in the quasiparticle continuum lying above 2A and 
heavily damped and partly because the important de- 
cay process of an amplitude mode into two H modes is 
beyond the reach of the random phase approximation. 

We shall pick up only the /i = v = component of 
the bubble diagrams involved in the vector interaction, 
i.e., Ily and n^, while discarding the longitudinal spa- 
tial components of the bubble diagrams whose cut-off 
dependence is not well controllable in the NJL model 
[see discussions around Eq. (|A3(I for details.] These spa- 
tial components, even if taken into account, would only 
change zero sound from superluminal to subluminal and 
would not essentially change the iJ-vector mixing and 
the behavior of attenuated zero sound in the CFL phase, 
as we shall see in the next section. 

We list the explicit expressions for H# , II^P , and IT^ in 
Appendix IA1 We shall omit the superscript "0" hereafter. 

It is convenient to arrange the polarizations and the 
coupling constants in a matrix form 



n 



-n M n v - 



G = 







G_ 

u 2 



(21) 



It is clear from this form that we can analyze the ff- 
vector mixing effect by picking up eigenmodes in the cor- 
responding channel. It should be noted that we always 
mean by the H channel the channel corresponding to 
phase fluctuations of the Cooper pair and that H as a 
physical excitation appears in an eigen-channcl given by 
a mixture of the H channel and the vector channel. 

The retarded propagator as sketched in Fig. [21 is ex- 
pressed in terms of II and G as 



1 



(2G)- 1 -n(w + 0+,g)' 



(22) 
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This propagator contains all information on the dynamics 
of collective excitations in the collisionless regime such as 
the mass, the attenuation width, and the energy disper- 
sion relation. As a simple exemplification, let us briefly 
make sure how the propagator (1221) is consistent with the 
presence of a massless H mode. In the limit of q = and 
lj — > 0, the mixed diagram II m vanishes and the denom- 
inator of the H propagator reduces to 6/G — 11^(0,0). 
Note that this denominator is proportional to the left 
side of the gap equation (JT^J. We thus find that it is 
zero given the explicit expression (| A2|) . This means that 
the propagator H22[l has a massless pole in the H channel 
corresponding to the Nambu-Goldstone boson. 

For the purpose of clarifying the overall properties of 
collective excitations, it is instructive to investigate the 
spectral function, which is defined as 

p(co,q) = --lmD R (uj,q). (23) 

In the subsequent sections we will calculate the dimen- 
sionless spectral functions in the H and vector channels: 

6 1 

P H= G HH> PV =G^ W VV> ( 24 ) 

where [p]hh is the H-H component of the 2x2 spectral 
matrix p and is the V-V component. We will then 

identify the poles of the propagator in the complex energy 
plane that characterize the spectral shape. 



IV. ZERO TEMPERATURE RESULTS 

In this section, we analyze the spectral functions at 
zero temperature. We clarify the properties of H phonons 
in the CFL phase and of zero sound in the normal and 
CFL phases. 



A. Phonon in a superfluid 

We start with an H phonon in the CFL phase. This 
mode, usually denoted as ff , ca n be interpreted as the 
Anderson-Bogoliubov mode |29j as encountered in neu- 
tral superfluids such as superfluid 3 He and superfluid 
neutron matter. In this mode, fluctuations in the phase of 
the superfluid order parameter are coupled with baryon 
density. Consequently, this mode behaves as a longitudi- 
nal compression mode; at zero temperature, its velocity 
is equal to the velocity of first sound. In the BCS frame- 
work, this mode can be viewed as a coherent superposi- 
tion of Cooper pair excitations with total momentum q, 
and hence it does not break up phase coherence in the 
superfluid state. We remark that as far as excitations as- 
sociated with baryon number are concerned, symmetric 
CFL quark matter behaves like a neutral superfluid even 
though each quark pair carries electric and color charge. 




FIG. 3: Spectral function in the H channel at zero tempera- 
ture. 

As we have already seen, the H propagator has a mass- 
less pole. Correspondingly, the spectral function natu- 
rally has a 5-function peak. In the H channel the mixing 
effect makes only a little difference in the spectral func- 
tion. As we shall see in the next subsection, however, it 
plays an important role in the vector channel. 

From Fig. |3| we can clearly observe that a massless 
H mode appears as a ^-function peak (strictly speaking 
there is a small width coming from a finite regulator for 
numerical computations) and that the continuum starts 
at the threshold lj = 2A above which individual pair ex- 
citations are kinematically allowed. The H mode always 
lies below the threshold; if lying above the threshold, it 
would be superseded by an excited Cooper pair, which 
eventually decay into two octet quark quasiparticles hav- 
ing an energy gap A. We remark that there appears a 
marked structure in the continuum around lj — 4A. This 
energy corresponds to a threshold above which individual 
pair excitations of singlet quarks are allowed to occur. 

We can also observe that two spectral peaks near 
u> = 2A and lj = 4A are appreciable for q > 120 MeV and 
that as the momentum q increases, they become sharper 
and sharper. This spectral enhancement of the contin- 
uum around lj — 2 A is analogous to the spectral enhance- 
ment as intensely argued in the chiral system of 7r and 
<7 mesons in connection with chiral restoration |30i l3lj | . 
In the present case no phase transition occurs with in- 
creasing but it seems reasonable to consider that a 
superfluid behaves more like a normal fluid for larger q, 
which will be illustrated in the next subsection. In other 
words, the medium in which H propagates undergoes 
continuous changes from superfluid type to normal fluid 
type with increasing q. 

This continuum structure characterized by the two dis- 
tinct gaps for octet and singlet quarks is a unique feature 
of the CFL phase, which can be seen neither in superfluid 
3 He nor in the 2SC phase, and provides useful informa- 
tion for identifying the CFL phase. 

We remark that the two peaks near lj = 2A and lj = 
4 A do not move substantially until q exceeds 2/i r . This 
is because the momenta of two quark quasiparticles, p 
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FIG. 4: The energy dispersion relation of H as a function of 
q at zero temperature. 



and \p + q\ , can sit on the Fermi surface simultaneously 
as long as q < 2/i r . We have confirmed that the peaks 
move upward for q > 2p r , although we will not present 
the results here. 

Another important property that can be observed from 
Fig. is the excitation energy of H as a function of q, 
i.e., the energy dispersion relation of H. The excitation 
energy of H is proportional to q at small q, and mono- 
tonically increases with q in such a way as to approach 
the threshold u) — 2 A asymptotically. This behavior is 
consistent with the QCD analysis performed by Zarembo 
in weak coupling |32| . For further clarity, we plot our nu- 
merical results for the energy dispersion relation in Fig.0] 
In the numerical calculations of the energy of H, we have 
ignored corrections by the mixing effect, which are in fact 
small at low energies and vanish at lo = 0. The dispersion 
relation shown in Fig. 0] thus fulfills 6/G— n#(cj, q) = 0. 

The earlier analyses based on the chiral effective La- 
grangian [33j are valid for u> and q much lower than 2 A. 
The present NJL model calculation is expected to repro- 
duce the results from the low energy effective theories. As 
can be seen from Fig. 21 the speed of H is vh = 1/ Vo, 
which is in agreement with the generally accepted value 
in the chiral effective Lagrangian approach. This is equal 
to the velocity of first sound, which is consistent with the 
fact that the restoring force of this mode is purely kine- 
matic. 

More technically, as shown in Ref. [3J|, the H velocity 
of l/\/3 originates from the dimensionality (three spatial 
directions and one temporal direction). In the rest of 
this subsection we shall calculate the H decay constant, 
Jh , by following the line of the standard NJL calculation 
of f„. Then, we shall consider how vh = l/\/3 results 
from the ratio between the decay constants in the spatial 
direction, f s H , and in the temporal direction, p H . 

We can obtain from the loop diagram shown in 
Fig. [S] as in the case of in chiral dynamics. The dia- 




FIG. 5: Diagram contributing to fu- 



gram can be expressed as 



«ftjff 



toT-j=FyS(po + qo,p+q)ig H qqT D S(po,p), 



(2tt) 4 2^6 



(25) 



where the factor 1/2 allows for the Nambu-Gor'kov 
counting, and the normalization factor 1/2V6 in front 
of Ty is determined so as to satisfy the proper current 
algebra. The iJ-quark coupling constant, gHqq, can be 
derived from the residue of the H propagator; the ex- 
pression for gHqq is listed in Appendix ^ Then, we 
find f H = (-1/a/6)A3^. [The sign of f H is irrele- 
vant to observables because onl y it s square, f H , appears 
in low energy effective theories [23].] After all we obtain 
the Goldberger-Treiman relation for quarks in the CFL 
phase, 



(fH ■ gHqq) 2 



(26) 



which was also derived in Ref. [34| from QCD in weak 
coupling. We note that gnqq remains finite and thus fn 
vanishes in the limit A — > as it should. 

When q^ is infinitesimal and the temperature is much 
lower than the quark chemical potential, it is sufficient 
to limit our deliberation here to the particle part of the 
quark propagators in Eq. 1251) . Then, one can show that 
the integrand of Eq. I|25|) is proportional to qo in the limit, 
q = and qo ~ 0, relevant for while this qo is replaced 
by {p-q)Pi in the limit, qo = and ft — 0, relevant for /J^ . 
By noting the three-dimensional rotational symmetry, 
one can rewrite J d 3 p (p-q)pi ... as (ft/3) J d 3 p .... Con- 
sequently, the ratio f H /f H amounts to g s Hqq /3g Hqq . By 
using this result, the relation v H — (f H / fn) 2 from the 
chiral effective Lagrangian, and the Goldberger-Treiman 
relation (|26|l . we finally obtain 



(27) 



From the derivation of this result, it is clear that 1/3 
comes from the spatial dimensionality. 

Let us next derive the expression for fn in the NJL 
model. Using the Goldberger-Treiman relation (|26|l and 
the approximate zero-temperature expression (|A8|) for 
g\j qq as derived in Appendix ^ we can express the tem- 
poral decay constant at zero temperature as 



t \2 



(Ph) 



Zn 2 r 
8tt 2 



(28) 
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FIG. 6: Spectral function in the vector channel, calculated in 
the normal phase at zero temperature. 



This result is the same as the one derived from QCD in 
weak coupling |33j |. Here we emphasize that the mean- 
field quark propagator in the NJL model takes the same 
form as that in QCD in weak coupling. As long as dia- 
grams composed of quarks are concerned, therefore, the 
NJL model and QCD yield essentially the same results 
within the mean-field approximation. This fact can ex- 
plain why the NJL model calculations of f v for CFL pi- 
ons have been found to be close to the QCD results in 
Refs. (3^, |3(j. We remark that the same kind of agree- 
ment between the mean-field NJL and QCD calculations 
can be also found in the parameters characterizing the 
Ginzburg-Landau free energy (24[. 

The results for fu and vh obtained in this subsection 
imply that the NJL model calculations should encompass 
the results from chiral effective Lagrangian approach and 
even the QCD results in the mean-field approximation. 



B. Zero sound 

We turn to the analysis of the spectral function in the 
vector channel at zero temperature. There are two im- 
portant features to be clarified in this subsection: the 
appearance of zero sound and the mixing effect between 
the H and vector channels. Zero sound is distinct from 
the H mode in the sense that it coherently involves a 
number of particlc-holc excitations rather than pair ex- 
citations. Wc first articulate how zero sound emerges in 
the vector channel in a normal fluid of quark matter. We 
then consider zero sound in the CFL phase where the 
mixing effect plays an important role. 



where the spectral function in the vector channel is plot- 
ted. In calculating the vector spectral function in the 
normal phase, we have fixed all the parameters at the 
values adopted in the previous section except A, which 
is set to be zero. We can observe from the figure that 
the continuum ranges < u> < q and that an undamped 
mode (zero sound) stands in the vicinity of the contin- 
uum, although it is rather hard to see in the plot that 
the (5-function peak is away from the continuum. 

The continuum in the spectral function stems from in- 
dividual excitations of quark particle-hole pairs. The 
continuum corresponds to the space-like region u) < q. 
In the present NJL model, the energy of a quasiparticle 
of momentum k is k as in the case of a massless free quark 
gas. For a real pair of a particle of momentum k\ and 
a hole of momentum £2, therefore, we can calculate the 
excitation energy as u> — (fci — /i r ) + (/i r — k?) < q, where 
q = ki — k\ is the momentum of the pair. 

If zero sound lies in the continuum, it would suffer Lan- 
dau damping. In this case, zero sound would be absorbed 
by a quasiparticle with momentum below the Fermi sur- 
face, which subsequently would be scattered into a state 
above the Fermi surface. This process corresponds to the 
decay of zero sound into a real particle-hole pair, which is 
kinematically allowed only when to < q. As we shall show 
explicitly, however, the calculated dispersion relation lies 
above the continuum. 

We now proceed to determine the velocity of zero 
sound by following a line of the conventional field- 
theoretical argument |37j . Let us focus on the expression 
for Ily, which is given by Eq. (|A3I) in the limit of A — > 0. 
At zero temperature the quasiparticle distribution func- 
tion vanishes. The dominant term that remains in IT^ is 
then 



n v (w,g) =S 18 



d 3 p 



(2tt) 3 

Q(Hr-p)Q(\p+q]- fj, r ) 



Q(p- fj, r )Q(fj, r -\f+ q\) 

\\p+q\-Hr\ + \p-pb 



{\\P+<i\- Vr\ + \P~ Vr\) 



(29) 



Here O denotes the Heviside's step function. The first 
(second) term in the square brackets corresponds to a 
virtual excitation of a particle (hole) with momentum p 
and a hole (particle) with momentum p + q. Because 
of these terms, the dominant contribution to Hy comes 
from a regime that satisfies p ~ \x r and u),q -C \i r . We 
can thus use the approximation, \p + q\ ~ p + p- q. Even- 
tually, can be estimated as 



1. In the normal phase 

In a normal fluid, where A = and the baryon num- 
ber is conserved, no mixing between the H and vector 
channels arises, and we can deal with the vector chan- 
nel alone. Our primary results are illustrated in Fig. 



IL v (oj,q) 



(2tt) 3 



dn- 



ln-q) 



Q (n ■ q) 



(30) 



where J dfl = 2tt fid cos 9 and n ■ q — cos 6. 

By defining the zero sound velocity as v = uj/q and 
using Eq. I|3U|I , we can rewrite the zero sound dispersion 
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relation, X/Gy + n^(a>, q) = 0, as 



wo 



In 



1 



- 1 = 



1 

Gy 



9/i 2 . 



(31) 



This agrees with the conventional form of the nonrela- 
tivistic zero sound equation (see Eq. (2.26) of Ref. [23). 
Since the right side of Eq. I|31[) is 2.304 for our parameter 
choice, we obtain vo = 1.0028 as a solution to Eq. (|31J) . 
Our numerical results for q = 40 MeV, on the other hand, 
has a S- function peak at ui = 40.11 MeV, from which we 
can derive the speed as v = 40.11/40 = 1.0028. This 
value is in excellent agreement with the analytic esti- 
mate. 

It is to be noted that the velocity of zero sound is 
slightly larger than unity. This uncausal behavior is due 
to our negligence of the longitudinal spatial components 
of the bubble diagrams in the vector channel, which con- 
tain a term proportional to the cut-off squared A 2 . We 
have confirmed that our calculations, if including the 
longitudinal spatial components, would give rise to zero 
sound of velocity below unity, which is damped in the 
continuum. However, it is more important to note that 
the Fermi velocity is unity since our model ignores the 
quark wave-function renormalization. Undamped zero 
sound, if present, would thus be automatically superlu- 
minal. In the high density limit of QCD in which one- 
gluon exchange dominates the quark-quark interactions, 
on the other hand, the bubble diagram in the vector 
channel yields a continuum in the spectral function at 
to < VFq < q, where vf is the Fermi velocity, and an 
undamped zero sound mode of velocity less than unity 
but larger than vf 0- Our calculations, therefore, are 
successful in reproducing the weak coupling behavior of 
zero sound predicted from QCD, except for the slightly 
superluminal velocity of sound. We thus believe that our 
framework gives a sufficient basis in analyzing the prop- 
erties of zero sound in the CFL phase, as will be discussed 
below. 




FIG. 7: Spectral function in the vector channel, calculated in 
the CFL phase at zero temperature. 
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FIG. 8: Spectral functions for q — 80 MeV in the vector chan- 
nel, calculated in the normal phase (dashed line), in the CFL 
phase without mixing (dotted line), and in the CFL phase 
with mixing (solid line) . The zero sound peak at u> ~ 80 MeV 
in the normal phase shifts to uo ~ 94 MeV with an attenuation 
width in the CFL phase. 



2. In the CFL phase 

Let us now examine zero sound and .ff-vector mixing 
in the CFL phase. We first present our numerical results 
for the spectral function in the vector channel and then 
discuss the underlying physics of the observed spectral 
suppression at small q and attenuation of zero sound. 

We plot the spectral function in the vector channel in 
Fig. [7| which is to be compared with Fig.^J The spectral 
functions displayed in Figs.EJandHare markedly different 
in two respects: In the CFL phase a (S-function peak 
corresponding to an H phonon in the space-like (u> < q) 
region appears, and a peak corresponding to zero sound, 
located atw~ 78 MeV for q = 60 MeV, at u ~ 94 MeV 
for q = 80 MeV, and at u ~ 111 MeV for q = 100 MeV, 
has a width in the continuum that extends above u = 2A. 
As a result of the .ff -vector mixing, the H phonon peak 
and continuum similar to the ones in the H channel as 



shown in Fig. |3| do appear in the vector channel, but 
have a vanishingly small spectral weight at q ~ 0. On 
the other hand, the observed width of the zero sound 
peak indicates the attenuation of zero sound. In order 
to see how the attenuation and the mixing effect affect 
the spectral function, we plot in Fig. |H1 three spectral 
functions calculated in the normal phase, in the CFL 
phase by ignoring the mixing effect, and in the CFL phase 
by including the mixing effect. 

First, we discuss the strong spectral suppression in the 
vector channel at q ~ and the underlying mixing effect 
between the H and vector channels. In the long wave- 
length limit q = 0, a clear distinction between particle 
and hole excitations is lost in the superfluid medium, 
and only the H (particle-particle) channel contributes to 
the spectral weight. For q far larger than A, on the other 
hand, the medium behaves more like a normal fluid rather 
than a superfluid, and the vector (particle-hole) channel 
exhibits a similar spectral function to that in the normal 
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phase. 

We can explicitly show the spectral suppression by not- 
ing that the bubble diagrams in the limit q — ► satisfy 
the relations, 

n ff (w,o) = riff (0,0) + — n v (w,o) 

= ^ + £^lv(^,o), (32) 

IL M (Lj ) 0) = -^a v {cj,0). (33) 

As a result of these relations, the vector component of 
the propagator matrix (|22[1 reduces to —Cry, which is ex- 
actly the tree- level value (the cross diagram in Fig. [5J. 
This indicates that the vector channel no longer contains 
the quark loops responsible for collective and individual 
excitations, leading to a vanishing spectral function. In 
other words, any excitations with small momenta in the 
vector channel tend to be absorbed into the H channel in 
the superfluid medium. As q becomes larger, the vector 
spectral function grows in the medium in which parti- 
cle and hole excitations are no longer indistinguishable, 
leading to a zero sound peak. 

This difference in the vector channel at small q between 
the normal and CFL phases could serve as a possible di- 
agnosis of superfluidity of quark matter if the spectral 
function in the vector channel could be measured in the 
lattice QCD simulation. Quark matter could be consid- 
ered to be in the normal phase when zero sound is found 
at small q, while in the superfluid phase when strong 
spectral suppression or weak mixture with an H phonon 
is observed. One could distinguish between zero sound 
and an H phonon from the fact that the H phonon ve- 
locity is smaller than the zero sound velocity by a factor 
of ~ l/\/3 except for Fermi-liquid corrections. In fact, 
zero sound has been found in the lattice simulation based 
on a QCD-like model If the simulation involves a 

superfluid phase, the vector channel measurement would 
provide a clear signature of superfluidity. 

We next discuss the attenuation of zero sound in the 
continuum. In order to extract detailed information on 
zero sound, we search for a pole of the vector propa- 
gator in the complex energy plane, which is responsible 
for a weak peak in the spectral function in the vector 
channel. For simplicity we ignore the -ff-vector mixing 
and the antiparticle contribution in locating the complex 
pole position. This kind of pole search has been per- 
formed in the context of chiral restoration in a hot and 
dense medium |3l| and of a phase transition from the 
normal to the 2SC phase Q. Appendix 151 is devoted to 
exhibiting how to perform the analytic continuation into 
the second Riemann sheet where the attenuated mode 
would be found. 

In Figs. OH and EH we plot the pole position, u, as a 
function of q. We can observe from Fig. [5] that the real 
part of to almost linearly increases with q once q becomes 
larger than ~ 2A. Since there is no quark excitation at 
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FIG. 9: The real part of the zero sound pole in the CFL phase 
at zero temperature. 
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FIG. 10: The pole position of zero sound in the CFL phase 
at zero temperature. The curve is slightly distorted when 
it passes through the threshold energy Kecu = 100 MeV for 
individual pair excitations of singlet quarks. 

low lu due to a nonzero pairing gap, zero sound, a coher- 
ent superposition of particle- hole pairs, can only emerge 
in the continuum where quasiparticle modes are allowed 
to excite. This explains why the energy of zero sound 
has a gap and stays slightly larger than 2A for q — > 0. 

Figure lrtjl illustrates how the pole moves in the complex 
energy plane with increasing q. The width of zero sound 
attenuation decreases with q. This attenuation occurs in 
the superfluid medium in which zero sound is absorbed 
by a Cooper pair, which subsequently breaks up into two 
quasiparticles. Consequently, the decrease in the width 
with increasing q is consistent with the tendency that the 
medium behaves more like a normal fluid for elementary 
excitations with larger q. Recall that zero sound in the 
normal phase is undamped in the present analysis. 

We conclude this section by summarizing the zero- 
temperature results for an H phonon and zero sound in 
the CFL phase. H phonons are only low-lying longitudi- 
nal excitations and survive as long as lu < 2A, while zero 
sound manifests itself as an attenuated mode for lu > 2 A. 
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The former is consistent with the earlier results from low 
energy effective theories. The latter is the point that we 
first clarified in the context of superfluid quark matter. 
The presence of the phase and zero sound modes in dif- 
ferent regimes of u> can also be seen in superfluid 3 He [l8j j. 
We also discovered a strong spectral suppression in the 
vector channel at small q. We remark in passing that all 
these properties of the CFL state are robust even if the 
longitudinal spatial components of the bubble diagrams 
in the vector channel are taken into account. 
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V. FINITE TEMPERATURE RESULTS 

In this section we examine thermal effects on collec- 
tive excitations in the collisionless limit. The collisions 
between quasiparticles, which are expected to play a 
role in modifying the properties of the collective modes 
at finite temperatures, will be addressed in our future 
study. At finite temperatures, as we shall see, the H 
modes are Landau damped by thermally excited quark 
quasiparticles, while the behavior of zero sound is es- 
sentially the same as that obtained at zero temperature. 
Throughout this section we limit ourselves to the case at 
T = 0.8T C = 14.52 MeV, which is high enough for us to 
perceive the temperature effect in the H channel, while 
not being in the immediate vicinity of the critical point. 
At this temperature, the pairing gap and the effective 
chemical potential can be calculated as A = 16.4 MeV 
and fi r — 435.3 MeV. The temperature dependence of 
fx r turns out to be tiny; even just above T c , /j, r is greater 
than the zero temperature value only by 2%. 

We first consider zero sound, which, at zero tempera- 
ture, appears as a low-lying mode in the normal phase 
and as an attenuated mode above 2A in the CFL phase, 
as we have seen in Sec. UVBl The behavior of zero sound 
at T = 0.8T C and at T = is almost indistinguishable. 
This is because T = 0.8T C = 14.52 MeV is low compared 
with the Fermi temperature corresponding to the quark 
chemical potential ~ 500 MeV. In the normal phase, un- 
less the temperature is considerably higher than T c , no 
drastic change occurs in the behavior of zero sound from 
that shown in Fig. |SJ 

In contrast to the case of zero sound, H modes in the 
CFL phase are damped only at nonzero temperatures. In 
the case in which thermally excited quark quasiparticles 
are present, the H modes suffer Landau damping. This 
process can be viewed as the upscattering of a quasipar- 
ticle by absorbing a H phonon, which is characterized by 
the nonvanishing factor [rip — rip'] in Eq. JUJl. The simi- 
lar damping process of the phase modes has been consid- 
ered in the context of neutral BCS superfluids [3q|. The 
Landau damping of the H modes can be observed from 
Fig. ^2 in which the spectral function in the H channel 
is plotted at T = 0.8T C . We find that the H modes have 
a finite decay width in contrast to the zero temperature 
case as shown in Fig. [2 The threshold for pair excitations 
is located at u> = 2A = 32.8 MeV. We also note that as a 



FIG. 11: Spectral function in the H channel, calculated in 
the CFL phase at T = 0.8T C . 
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FIG. 12: The pole position of the diffusive mode responsible 
for the Landau damping in the CFL phase at T — 0.8T C . 



result of the Landau damping, the spectral function has 
a continuum below the H peak. 

Detailed analysis shows that the velocity of the H 
mode at T = 0.8T C is vh = 0.554 which is slightly smaller 
than This is because the superfluid part rather 

than the whole system participates in the phase fluctua- 
tions. 

Now, we briefly discuss the possible relevance of our 
results to the transport properties in the CFL phase. In 
the H channel, by setting II m = for simplicity, we can 
find a pole on the imaginary axis in the complex energy 
plane. This pure imaginary pole corresponds to a diffu- 
sive mode that does not propagate at all. We identify 
this imaginary pole as responsible for the low-lying con- 
tinuum due to the Landau damping. 

In Fig. El we plot the pole position as a function of 
the momentum q. The result can be parametrized as 
Imw = -30.9 - 5.11xl0~V [MeV]. This result suggests 
that H modes in CFL quark matter suffer diffusion with 
the diffusion constant 5.11 x 10~ 3 MeV -1 in addition to 
absorption. This diffusion coefficient arises even in the 
collisionless limit, but has to be modified by collisions 
between quasiparticles, which are more relevant at higher 
temperatures. 

Although description of collective modes in the hydro- 



13 



dynamic regime is beyond the scope of the present paper, 
we shall comment on what is expected in the hydrody- 
namic regime. In the normal phase, as in a normal liq- 
uid 3 He, zero sound in the collisionless regime should 
continuously transform to first sound in the hydrody- 
namic regime. In the superfluid medium, as discussed in 
Ref. [3^ and also pointed out in Ref . |2(j , the mode corre- 
sponding to phase oscillations in the collisionless regime 
should change into second sound in the hydrodynamic 
regime. It would be an interesting future problem to de- 
scribe such continuous transitions in the nature of sound 
modes within the framework of kinetic theory. 

VI. CONCLUDING REMARKS 

In this paper we have studied the longitudinal spec- 
tral functions in the H and vector channels in a super- 
fluid of CFL quark matter. We have elucidated the zero- 
temperature properties of the vector channel: the strong 
spectral suppression at small momenta due to the H- 
vector mixing and the zero sound attenuation above the 
energy threshold for pair excitations. The important im- 
plication of these results is that the spectral function in 
the vector channel is sensitive to the state of quark mat- 
ter. If lattice QCD simulation at finite density is feasible 
in the future, the vector channel measurements could be 
of great use for identifying superfluidity of quark matter. 

Extension of the present analysis to the case of nonzero 
quark masses is desirable for the description of elemen- 
tary excitations in realistic situations. In this case, as 
discussed in Sec. |n] the possible coupling of phase and 
baryon density fluctuations with electric and/or color 
charge density fluctuations could change the properties 
of the longitudinal collective excitations drastically. In 
the presence of gapless quark modes as in the gapless 
CFL phase [25| , the continuum of the quark excitations 
and hence the attenuation of the collective modes would 
be enhanced. 

In the present analysis we have utilized the NJL model 
for quark matter. The NJL model description turns out 
to reproduce the results for H phonons known from low- 
energy effective theories and even from QCD in the mean- 



field approximation. In the collisionless regime, however, 
the longitudinal spectral functions remain to be clarified 
from QCD even in weak coupling. Furthermore, collec- 
tive modes in the CFL phase have yet to be examined in 
the hydrodynamic regime. It would be instructive to de- 
scribe low-lying excitations of CFL quark matter in both 
regimes by using a relativistic version of the Landau the- 
ory of Fermi liquids |19|, which effectively encompasses 
the NJL model. 
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APPENDIX A: EVALUATION OF THE BUBBLE 
DIAGRAMS 

In this appendix, we present the explicit expressions 
for the polarization functions in the H , vector, and mixed 
channels. We also provide the calculation of guqq which 
is a derivative of the polarization with respect to either 
energy or momentum. Since the quark propagator is di- 
vided into the particle part and the antiparticle part as 
displayed in Eq. J7J) , the bubble diagram in the / channel 
(/ = H, V, M) can be divided into four parts as 

n ; = n pp + nf a + n pa + n ap , ( ai ) 

where the superscript "p" and "a" denote the particle 
and antiparticle parts of the quark propagator, respec- 
tively. 

We begin with the expression for the polarization Tin 
in the H channel. The particle-particle part can be writ- 
ten as 



II£V, q) = 2 / A ** [1 + P ■ (p+fl)] -^7 

J ) e A e A 



e A e A +(P~ Mr)(|P +q\ ~ Mr) + A 5 



£a £ a' - (P - Mr)(|P + q\ ~ Mr) - A 5 



uJ + e A + e A iv-e A -e A 



[l - n F - rip'] 



(A2) 



u> — e A + e A oj + e A - e 



A C A 



— |_ p "F 



+ - (same with A -> 2 A) . 
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where we have used the notations for the quark excitation energies, e A = e A (p) and e A ' = e A (p + q), and for the Dirac- 
Fermi distribution functions, rtp = [e e ^/ T + l] _1 and n^' = [e e A l T + The antiparticlc-antiparticle part, Il^f, 

can be readily obtained from expression (|A2fl by replacing fi r by — /j, r . [Note that by this replacement e~ is replaced 
by e + .] One can obtain the expression for the particle-antiparticle part 11^ by replacing, in Eq. (|A2|) . (p — /i,.) by 

— (p + /i- r ), e A by e A' n F n F' anc ^ (p + by 1 — p- (p + q). The antiparticle-particle part, n^ 5 , likewise 

obtains. 

The particle-particle part of the bubble diagram in the vector channel reads 



[l+p- (p + q)]- 



1 



d 3 p 

w c a c a 

-(p-Hr)(W+q\ -Mr) + A 2 



£ A £ a' + (P - fJr)(W+ Q\ - Mr) 



(A3) 



+ - (same with A -> 2A). 



This is the same as except for the sign in front of 
the factor (jp — n r )(\p + q\ — ^r) and the coefficient af- 
fixed to the singlet contribution. This change in the sign 
can be intuitively understood once the limit of A — > is 
taken. In this limit, as discussed in Sec. IIVBI the inte- 
grand is nonzero only when one of the factors (p — (i r ) 
and (\p + q\ — [i r ) is positive and the other is negative. 
This situation corresponds to a simultaneous excitation 
of a quark particle above the Fermi surface and a quark 
hole below the Fermi surface, rather than pair excitations 
which play a role in H^f ■ In the presence of the energy 
gap A, the clear distinction between particles and holes 
is lost since the quasiparticle state appears as a linear 
superposition of a particle state and a hole state. The 
other parts, Il^f, IIy a , and II^P, obtain in the same way 
as in the H channel. 

We remark that IIy P is almost independent of the 
cut-off A because the momentum integration around the 
Fermi surface contributes dominantly to IIy P . For suffi- 



ciently large p to satisfy |p+<7| > ^t r , a hole is not allowed 
to excite at zero temperature, so that the momentum in- 
tegration has a natural cut-off. In contrast, Ily 3, has an 
appreciable dependence on A since an antiparticle can 
arise at arbitrary p. This cut-off dependent contribution 
is, however, suppressed by a factor of q 2 /p 2 because of the 
orthogonality between the positive and negative energy 
states, i.e., trAp +q AJ = 1 — p ■ (p + q). If we take account 
of the spatial component of the polarization, Ily 5 , for 
instance, a particular choice of the polarizing direction 
leads to breaking of the orthogonality. As a result the 
particle-antiparticle contribution of Ily has a term pro- 
portional to the cut-off squared A 2 , and hence the cut-off 
dependence of Ily is not well controllable in contrast to 
the case of 11^ considered here. 

We turn to the expression for the polarization Hm in 
the mixed channel. The particle-particle part reads 




and , n^, and n^ p obtain in the same way as in the as 
H channel. 



From the expression for Hh obtained above, we can 
calculate the coupling constant between H and quarks, 
which is determined by the residue of the H propagator where m# is the H mass. We note that since mn = 0, 



an 



H 



dq^qu 



(A5) 
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the mixed part Hm does not contribute to gnqq- Since 
Lorentz symmetry is broken in a many particle system, 
the derivatives in the temporal and spatial directions lead 
to different results, which in turn yield the velocity of the 
H mode smaller than unity. 

We can write down the derivative in the temporal di- 
rection as 



dngV,o) 



=2 



Hqq) doJ 2 

A d 3 p l-2n F 



U1=0 



(A6) 



(2?r) 3 ( 3 +2(samewithA 



2A) 



for the particle-particle part. Here we note that 
IF^(ui, 0) = and that is suppressed compared with 
II^ 3 since e A in il^ is replaced by e A in 11^. At zero 
temperature we can estimate the integral in Eq. i|A6(l . by 
assuming that it is dominated by momenta close to the 
Fermi surface, as 



2 /-oo 



dp 



d 3 P 1 

(2tt)3 Je^) 3 ~ 2tt 2 (p 2 + A 2 ) 3 / 2 ~ tt 2 A 2 ' 

We thus obtain the temporal iJ-quark coupling constant 
at zero temperature as 



(9 Hqq) 



9tf 
4tt 2 A 2 



(A8) 



We next estimate from Eq. (|A5|) the spatial iJ-quark 
coupling constant at zero temperature. The differentia- 
tion in the spatial direction involved in this estimate is 
more complicated than that in the temporal direction. 
By expanding the integrand of 11^(0, q) with respect to 
q, we find that the term proportional to q 2 is smaller than 
the term proportional to (p ■ q) 2 by a factor of A 2 /fj, 2 . 
We thus neglect the former term. Then, the integrand of 
<9nPP(0, q)/dq 2 reduces to 

A 4 + (e A ) 2 ( M 2 - 3^.p + 3p 2 ) + A 2 (/i 2 - 2^ rP - 2p 2 ) 

(A9) 

Note that the integral with respect to p is dominated 
by momenta close to the Fermi surface. We can thus 
simplify the integrand by dropping the A term and by 
setting p = [i r in the numerator as 



^ 2 [( £A ) 2 -3A 2 ] 



(A10) 



We finally obtain 



(.9 



\-2 



dWg(0,q) 



Hqq) Q q 2 

A d 3 p (p-q) 2 _ 
(2ir) 3 ( £a )3 - 



g=0 



(All) 



^{9 Hqq 



a) - 



where we have used the relation J dp (e 
(2/3A 2 ) J e?p(e A )~ 3 , and the factor 1/3 originates from 
the average over three spatial directions. 











FIG. 13: Poles off the real axis in the complex p plane. We 
can perform the analytic continuation by picking up the pos- 
itive residue of poles in the upper-half plane and the negative 
residue of poles in the lower-half plane. 



APPENDIX B: ANALYTIC CONTINUATION TO 
THE SECOND RIEMANN SHEET 

Here we explain how to perform the analytic contin- 
uation to the second Riemann sheet on which possible 
unstable states would have a pole. In the imaginary-time 
formalism of finite temperature field theory, the real-time 
retarded propagator, D R (Lu,q), can be obtained by the 
analytic continuation, D R {uo,q) — D(qo = to + i0 + ,q). 
When u> is larger than the threshold for single particle 
excitations, the propagator has a finite imaginary part, 
which leads to the spectral function having a continuum. 
Then, the propagator and hence the polarization II (to, q) 
has a branch cut in the complex energy plane. 

If we take II(z = £ — irj, q) as an analytically continued 
function, this polarization in the limit of 77 — > is not 
continuously connected to U(( + i0 + , q) in the physical 
sheet due to a branch cut in the continuum region. In 
order to remedy this situation, it is convenient to add the 
discontinuity, 

DiscII(w, q) = U(oj + i0+, q) - U(uj + i0~,q) . (Bl) 

Once we know the analytic continuation of Discn(oj, q) 1 
we can locate the pole in the second Riemann sheet by 
using 11(2:, q)+DiscII(z, q) as a properly analytically con- 
tinued function. 

In our calculation II is not written as a closed expres- 
sion, but, as listed in Appendix IA1 an expression involv- 
ing the momentum integration with respect to p. For 
the energy lo + i0 + , the integrand has in general several 
poles stemming from 1/(lo — e A — e A ) and other simi- 
lar terms in the complex p plane. These poles lie either 
above or below the real axis as shown in Fig. 1131 Then, 
as is obvious from the figure, DiscII is nothing but the 
contour integrals around the poles. If a pole is above 
(below) the real axis, the contour goes in the positive 
(negative) orientation. Straightforward generalization of 
this computational rule to the case of complex to yields 
the analytic continuation. For complex uj — z, poles are 
not necessarily distributed in the vicinity of the real axis. 
Then, we can calculate DiscII(z,q) from the contour in- 
tegrals around the poles as shown in the lower part of 
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